Off-Equilibrium Effective Temperature in Monatomic Lennard-Jones Glass 
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The off-equilibrium dynamics of a monatomic Lennard-Jones glass is investigated after sudden 
isothermal density jumps {crunch) from well equilibrated liquid configurations towards the glassy 
state. The generalized fiuctuation-dissipation relation has been studied and the temperature depen- 
dence of the violation factor m is found in agreement with the one step replica symmetry breaking 
scenario, i. e. at low temperature m{T) is found proportional to T up to an off-equilibrium effective 
temperature %jf, where m{%jf) — l. We report %jf as a function of the density and compare it with 
the glass transition temperatures Tg as determined by equilibrium calculations. 



PACS Numbers : 61.20.Lc, 64.70.Pf, 02.70.Ns 

In the last few years the off-equilibrium dynamics of 
glassy systems has been the object of intensive analyti- 
cal and computational studies, due to the very rich phe- 
nomenology exhibited and the possibility to obtain de- 
tailed information on the phase-space properties of the 
glassy system itself. One of the most important fea- 
tures of the off-equilibrium regime is the generalization 
of the fluctuation-dissipation relation . At equilibrium 
the Fluctuation-Dissipation Theorem (FDT) ||] states 
that the susceptibility, XAB{t), that describes the re- 
sponse of the variable A{t) to an external field conju- 
gate to the variable B{t) is proportional, through the 
factor P={KbT)~^ , to the derivative of the cross cor- 
relation function of the variables A{t) and B{t) them- 
selves, CAsit). In the off-equilibrium case both suscepti- 
bilities and correlation functions are two-times quantities 
and the FDT no longer holds. However, as suggested by 
Cugliandolo and Kurchan through analytical studies of 
soluble spin glass models the FDT can be generalized: 
the ratio between XAB(ti,t2) and f3CAB{ti,t2) defines 
a function, X(ti,t2), that depends on the times ti,t2 
only through the actual value of the correlation function 
CAB{ti,t2), X(ti,t2)=X(CAB{ti,t2)). Thc bchavior of 
X{Cab) gives information on the equilibrium order pa- 
rameter in the low temperature phase , in agreement 
with numerical experiments on spin glasses Q . 

The conjecture of the similarity between structural 
glasses and some spin glass model has been the ba- 
sis to extend these analytical results on spin glasses to 
structural glasses. According to the latter conjecture 
one expects, in structural glasses, that X{Cab) is a 
two values function, with X{Cab)~^ at short times and 
X{CAB)^Tn < 1 in the non trivial long time region (the 
so called aging region) . The constant m is expected to be 
proportional to the temperature of the system up to an 
effective temperature Teff, above which m=l and, there- 
fore, the FDT is recovered. In the numerical tests of 
this scenario performed so far the system is brought 
off-equilibrium at constant density by sudden tempera- 
ture jumps [quench), and the effective temperature T^g 



is found to coincide with the glass transition temper- 
ature Tg. This is explained assuming that during the 
quench the system "captures" the properties of the 3A^- 
dimensional potential energy surface at that temperature 
where the relaxation time becomes comparable with the 
simulation time [Tg] and, at lower T, it is no longer able 
to equilibrate. In this way the actual aging dynamics co- 
incides with those at Tg even if the "true" temperature is 
T<Tg. This is a reasonable scenario, but it is still unclear 
what happens if the system is brought out of equilibrium 
without crossing the thermodynamic state [Tg , p) . 

In this Letter we numerically investigate the off- 
equilibrium fiuctuation-dissipation relation in a sim- 
ple monoatomic glass, whose particles interacts via a 
Lennard-Jones (LJ) potential modified in such a way 
to avoid the crystalization occurring in standard LJ. 
At variance with the usual constant density tempera- 
ture jumps method (quench), we use density jumps at 
constant temperature (crunch) in order to produce off- 
equilibrium configurations. As expected, at low temper- 
ature, we find violation of the equilibrium FDT in agree- 
ment with previous computational results in other model 
glasses ^] . We also find a linear temperature dependence 
of the off-equilibrium parameter m in the aging region, 
thus confirming the conjecture of one step replica sym- 
metry breaking scenario for real glasses [Q . More impor- 
tant, the effective temperature T^jf, below which m<l, is 
found very close to the glass transition temperature Tg 
at the density reached after the jump. While in the usual 
quench from high T at constant p it seems natural to vi- 
sualize the process of freezing at Tg during the quench, in 
our crunch method the identification of the two tempera- 
ture Tejf and Tg evidences how the effective temperature 
depends only on the density of the final state and does 
not depends neither on thc initial state nor on the path 
in (T, p) plane. This observation gives a strong support 
to the idea that the aging process can be considered as 
a pathway in the 3A^-d energy landscape from the off- 
equilibrium situation towards the glassy minima . 

The investigated system is a Modified Lennard-Jones 
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(MLJ) model glass, with N=256 particles. The poten- 
tial energy is V=Vlj+SV, where Vlj is the usual 6-12 
Lennard-Jones interaction (in the following all the di- 
mensional quantities are expressed in reduced units) and 
dV is a many-body term that inhibits crystalization: 



6V = aS,- 9{S{q) - So) [S{q) - So] 



(1) 



Here S{q) is the static structure factor and the parame- 
ters in SV (a*=0.8 and 5*0=10) have been tuned in or- 
der to avoid crystalization without introducing a per- 
turbation on the "true" LJ dynamics ||T^. The sum 
is made over all q with q„iax—^<\<l\<<lmax+^, where 
q^^^=7.12{p*y/^ and A*=0.34. With this potential we 
can investigate all the phase space without falling down 
in spatially ordered configurations and with negligible 
corrections to the equation of state of monatomic LJ. 

First we study equilibrium properties of the system, 
determining the liquid-glass transition line in the (T, p) 
plane. The system is equilibrated at high temperature at 
different densities, from p*=0.87 to p* = 1.21, then it is 
slowly cooled {dT* /dt*=4:.2-10~'^) measuring the poten- 
tial energy as a function of T. 




FIG. 1. Potential energy vs temperature (all in LJ units) 
during the cooling process from high temperature. The open 
circles are the molecular dynamics data, the solid line is the 
fitted curve (see text) and the dashed line is the harmonic 
extrapolation at high temperatures. 

In Fig. |l| we report an example of the potential energy 
vs T* at density p*=0.99. It is evident the existence of 
a temperature T* that marks the transition between two 
smooth regimes: at high T, in the liquid, as conjecture 
by Rosenfeld and Tarazona the potential energy fol- 
lows a T^/^ low while at low T, in the glass, the linear 
behaviour expected for harmonic glasses is recovered. In 
this way we are able to reconstruct the liquid-glass tran- 
sition line T*{p*), as reported in Fig. ||. Below p*=0.87 
a spinodal decomposition takes place, evidenced by the 
appearance of "bubbles" in the sample. 



We focus now on the fluctuation-dissipation relation. 
This relation connects the correlation function to the 
response to an external field. Let be A{t) and B{t) 
two generic microscopic quantities, and lS.TL=\h{t)B{t) 
a perturbation added to the Hamiltonian, with h{t)—0 
for t<tw and A an adimensional control parameter. If 
we define the correlation function at zero perturbation, 
CAB{t,tyj), as: 

CAB{t,t^)^{A{t)B{t^))^^,, , (2) 

and the corresponding response function xab- 



XAB\t,tw) = hm -- 



(3) 



(hereafter the subscript A indicates that the aver- 
age is performed in presence of the perturbing term 
in the Hamiltonian) then the equilibrium fluctuation- 
dissipation relation takes the form: 



XAB{t,tw) = dCABit,tw)/dt^ 



(4) 



We notice that even if at equilibrium the presence of 
is redundant, as XAB{t,tw) and CAB{t,tw) depend only 
on t — tyj, we choose to use a two times formalism in order 
to easily generalize the above formulas to off-equilibrium 
case. 
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FIG. 2. The liquid-glass transition line in the plane (T* , p*) 
in reduced units. The black points are determined from 
crossover in the temperature dependence of equilibrium po- 
tential energy (the dashed line is a eyes guide). The two grey 
squares are off-equilibrium effective temperatures Tejf deter- 
mined from fiuctuation-dissipation relation. Arrows qualita- 
tively describe the crunch process in the off-equilibrium anal- 
ysis (from Pq—0.95 to pl=1.14 and to p2 = l-24). 

Introducing the integrated response function 
Rab (t, tw), 



RAB{t,tw) = Jt dt' XAB{t,t') , 



(5) 



we obtain from eq. (^) the fluctuation-dissipation relation 
in the integrated form: 
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RAB{t,t^) = l3[CAB{t,t^) - CAB{t,t)] . (6) 

In the present work we use the quantities A{t) and B{t): 

A{t) = V2 (gA^)-iE. sin(<7l • nit) + ct>i) , (7) 
B{t) = -V2 g-iS, sin(g' • n{t) + 0,) , (8) 

where is a uniform random phase, qi = qsi, s,; a vector 
whose Cartesian components are random variables ±1, 
and (7=27rn/L, with n an integer and L the sample size. 
Moreover we choose h{t) = hoO{t — tw) {hl = l hereafter). 
From Eq. (H), after an average over the random phases 
(pi, one gets: 

CABit,t^) - -iq^N)-^^, { COs(q, • Sn{t,t^)) , 

(9) 

where 8ri{t,tw)=ri{t)-ri{tw), and, according to Eqs. (H) 
and d), 

RAB(t,t^)^V2{q\N)-^T.,[ ( s\ii{q, ■ n{t) + ), + 

- (sin(gl.r;(<) + 0,)>..o] • (10) 

(Fortunately in the limit N ^ oo the r.h.s of the equa- 
tions becomes self-averaging). A more suitable and easy 
to compute form is obtained in the limit qi ■ Svi In 
this limit, and averaging over s^, Eq. (j^) becomes: 

CAs(t,t^)^-9"' + A(t,i^)/2 , (11) 

having defined the mean square displacement: 

A(i,t,,) =iV-iS, (|r«(i) -^^(t-.)P>.=o- (12) 
Similarly, Eq. ( [lo| ) becomes: 

RAB{t,t^)c^{XNr^^,[{h-n{t)),~{f,-n{t)),^^] (13) 

where we have introduced the force A/i acting on the par- 
ticle i due to the perturbation term in the Hamiltonian: 



R{Cab) = -P J, 



CAB{t,t) 



dCAB X{Cab) , (17) 



fi = V2si cos(gi • fi{tw) + 4>i) 



(14) 



this is a random variable due to the randomness of (f, 
and Si. Therefore, from Eqs. ([ll| , [l3| ), Eq. (^ becomes 



2 R{t,t^) ^ 13 A{t,t^) 



(15) 



In the off-equilibrium case one can always define a vio- 
lation factor X(t, tw) in such a way to rewrite the Eq. (^) 
as follows: 

XAB{t,t^) = -(3 XAB{t,t^) dCAB{t,t^)/dtyj . (16) 

In some recent papers [|| it has been conjectured that, in 
the aging region, XAB{t,tw) depends on its arguments 
only through the correlation function CAB{t,tw)- This 
allows to generalize the fluctuation dissipation relation 
(i) as: 



dR{CAB)/dCAB = P X{Cab) ■ 
In our specific case we obtain: 

2 dR{A)/dA = 13 X{A) , 



(18) 



(19) 



which generalizes Eq. ( |l5|) . In this work we have studied 
the generalized FDT, in the form expressed by Eq. (|is|). 
For sake of semplicity rather than perturbing the system 
with the random force fi in Eq. (14), we have used a 
random force fl with the same variance of fi ((/^) — 
1) and Cartesian components ±1 with equal probability. 
Equation ( p^ ) has been numerically tested at equilibrium 
(see the upper part of Fig. 

The off-equilibrium states are obtained by crunches 
from an initial equilibrated configuration at density 
Po=0.95 at a given temperature T to a final state at the 
same T and higher density (p*=1.14 and p2=l-24). We 
set the time of the crunch as t=0. For each investigated 
T and p we perform an isothermal molecular dynamics 
simulations with A=0 and calculate A{t, t^) as a function 
of t and Starting from the same t^O configurations, 
after a waiting time we switch on the external field 
(A^^O) and measure the response R{t,tw). The strength 
of the external field A*=0.27 was chosen after extensive 
tests of the linear dependence of the response on the per- 
turbation. All the measured quantities are averaged over 
50 initial configurations and 10 random extractions of 
the variables The values investigated are t*^ = l, 
5, 10. As an example, in Fig. ^ we show the quantities 
A and 2R/ (3 after a crunch from pg=0.95 to p2~l-2A for 
three different temperatures, r*=1.76, 0.96, 0.48, and 
for t^=5. In the left hand part the two quantities are 
plotted versus log(i/t^) while in the right side 2R/P is 
plotted as a function of A. All the data are compatible 
with the assumption that violation factor X depends on 
t and only through A(t, t^), as verified comparing the 
results obtained for different values of the waiting time. 
Moreover the measured values AT (A) are well represented 
by a piecewise constant behavior (A=l for short times 
and A<1 at longer times) suggesting a one step replica 
symmetry breaking scenario for structural glasses as con- 
jectured in and already found in |^ in the cases of 
quenching systems of soft spheres and LJ binary mix- 
ture. 

In Fig. ^ we show the T dependence of the violation 
factor m{T) in the aging region, for the two final densities 
analyzed. It is evident a linear behavior up to a tempera- 
ture Tejf , as found at low temperature in all know cases in 
which one step replica symmetry holds. Fitting the data 
with the expression m{T)^aTelTejf - T)+0{T - T^g) we 
extract the corresponding Tf,ff values, which turn out to 
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correspond with the previously measured equilibrium val- 
ues of Tg as evidenced Fig. |l|. The coincidence between 
Tfiff at a given final density and Tg at the same equilib- 
rium density evidences supports the hypothesis that, in 
the aging region, the point representing the system in the 
3-/V-d energy landscape moves from the region pertaining 
to "high temperatures", where it has been brought by 
the sudden density change, towards those configurations 
where the characteristic time for structural rearrange- 
ment "diverges" P|. 
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FIG. 3. Mean square displacement A and integrated re- 
sponse function 2_R//3 in reduced units at p2 = 1.24 for three 
temperatures r*=0.48, 0.96, 1.76. The left side shows the log 
time dependence of the two quantities. In the right side 
the response function vs A. Dashed lines indicate equilib- 
rium fluctuation-dissipation relation, while full lines fit the 
ofF-equilibrium aging region. 

In conclusion we have, for the first time, studied the 
off equilibrium dynamics of a simple monatomic glassy 
system: a LJ potential where the crystalization has been 
inhibited adding a small many-body term. The system 
has been brought off-equilibrium via isothermal density 
jumps and the generalized fluctuation-dissipation rela- 
tion has been studied. The found piecewise constant 
behavior of the the violation factor and its linear tem- 
perature dependence agrees with one step replica sym- 
metry scenario conjectured for structural glasses. For 
a given density, the emergent off-equilibrium effective 
temperature T^-g seems to be very close the glass tran- 
sition temperature Tg , determined by equilibrium molec- 
ular dynamics measures. Then, even if brought out of 
equilibrium by crunches at constant T at final density p 
without crossing the Tg{p) transition point, the system 
captures the properties of energy surface at T^jf ^ Tg [p) . 



The off-equilibrium effective temperature T^g is then in- 
dependent on the particular path in (T, p) plane crossing 
the transition line, and depends only on the final reached 
density. 




FIG. 4. Temperature dependence of the violation factor 
in the off-equilibrium aging region for two densities: pi=1.14 
(open circles) and p2 = 1.24 (full squares). From the fits 
(dashed lines) one obtains: Xlj (pi)=0.86 and IJj (p2) = l-43. 
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